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By cyclic adiabatic change of two control parameters of an optical trap one can induce a circulating 
current of condensed bosons. The amount of particles that are transported per period depends on 
the "radius" of the cycle, and this dependence can be utilized in order to probe the interatomic 
interactions. For strong repulsive interaction the current can be regarded as arising from a sequence 
of Landau-Zener crossings. For weaker interaction one observes either gradual or coherent mega 
crossings, while for attractive interaction the particles are glued together and behave like a classical 
ball. For the analysis we use the Kubo approach to quantum pumping with the associated Dirac 
monopoles picture of parameter space. 



I. INTRODUCTION 



Dimers and Trimers 



Understanding the complicated behavior of quantum 
many-body systems of interacting bosons has been a ma- 
jor challenge for leading research groups over the last few 
years. In fact, the growing theoretical interest was fur- 
ther enhanced by recent experimental achievements. The 
most fascinating of these was the realization of Bosc- 
Einstein condensates (BEC) of ultra-cold atoms in op- 
tical lattices (OL). This allows for the vision that the 
emerging field of atomtronics (the atomic analog of quan- 
tum electronics) will result in the creation of a new gen- 
eration of nanoscale devices. A major advantage of BEC 
based devices, as compared to conventional solid-state 
structures, lies in the extraordinary degree of precision 
and control that is available, regarding not only the con- 
fining potential, but also the strength of the interaction 
between the particles, their preparation, and the mea- 
surement of the atomic cloud. The realization of atom 
chips EL "conveyor belts" [|J , atom diodes and transis- 
tors [3, |j| is considered a major breakthrough with po- 
tential applications in the field of quantum information 
processing Q , atom interferometry [f| and lasers [l], E} ■ 



The theoretical and experimental study of driven dy- 
namics in a few site system using optical lattice tech- 
nology [l3[ is state of the art. So far mainly single and 
double site (dimer) systems were in the focus of actual re- 
search. The study of a three-site (trimer) system adds an 
exciting topological aspect which we would like to explore 
in this work: the possibility to generate in a controlled 
way circulating atomic currents. 

The possibility to induce DC currents by periodic (AC) 
modulation of a potential is familiar from the context of 
electronic devices. If an open geometry is concerned, it is 
referred to as "quantum pumping" |14| . while for closed 
geometries we use the term "quantum stirring" [l5j]. In 
the present paper we consider stirring of condensed par- 
ticles [HI in a few-site system which is described by the 
Bose-Hubbard Hamiltonian (BHH). 



Recently, BECs in driven optical lattices have received 
a lot of attention. It was pointed out Q that the study of 
the energy absorption rate (EAR) can be used to probe 
the quantum phase of the BEC. Specifically, the excita- 
tion spectrum of the system was determined by measur- 
ing the EAR induced by a periodic modulation of the 
lattice height. The position and height of the peaks of 
the EAR versus the driving frequency give valuable in- 
formation about the interaction strength and incommen- 
surability of the system. The experimental activity on 
EAR spectroscopy and its applications triggered a theo- 
retical interest in understanding and predicting the EAR 
peaks @. Simultaneously, big efforts were dedicated to 
the study of driven dynamics of smaller optical lattices 
like single and double site systems [J, Ell EH EH aiming 
to understand how to tame quantum dynamics. 



B. Main observation 



The controlled stirring operation produces an adia- 
batic DC current in response to a cyclic change of con- 
trol parameters of the optical potential that confines the 
atoms. We find that the nature of the transport process 
depends crucially on the sign and on the strength of the 
interatomic interactions. We can distinguish four regimes 
of dynamical behavior. For strong repulsive interaction 
the particles are transported one-by-one, which we call 
sequential crossing. For weaker repulsive interaction we 
observe either gradual crossing or coherent mega cross- 
ing. Finally, for strong attractive interaction the particles 
are glued together and behave like a huge classical ball 
that rolls from trap to trap. 
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FIG. 1: (Color online) Illustration of the model system. Ini- 
tially, all particles are located on the upper site (i = 0) which 
represents the "shuttle" . In the first half of the cycle (a) the 
on-site potential vo = e is raised adiabatically slow from a 
very negative initial value and the particles are mainly trans- 
ported via the fci bond to the "canal" which is represented 
by the strongly coupled canal sites (i = 1,2). In the second 
half of the cycle (b) the bias in the coupling is inverted and 
the particles are mainly transported back from the canal to 
the shuttle via the bond. 



C. Scope and outline 

The present paper has several objectives: (i) To intro- 
duce an illuminating picture for the analysis of transport 
in a few-site system based on the adiabatic formalism, 
(ii) Using this picture to argue that there are four dis- 
tinct dynamical regimes dependent on the strength of the 
interactions, (iii) To propose a controlled stirring pro- 
cess that can be utilized in order to probe the strength 
of the interactions in a trimer system. 

The paper is organized as follows: The quantum trimer 
is introduced in Section 2. We then describe qualita- 
tively the stirring process (Section 3), and briefly review 
the Kubo formula approach to quantum pumping fl7j 
(Section 4) which is based on the theory of adiabatic pro- 
cesses [l8| (see Appendices A,B). In Section 5 we describe 
a reduction of the three-site Hamiltonian that allows us 
to regard the stirring of one particle in a trimer as a 
Landau-Zener (LZ) crossing in a two-level system. The 
analysis is extended to many particles in Section 6 where 
the variation of the adiabatic energy levels as a func- 
tion of the control parameters is analyzed. We regard 
the transport as a sequence of LZ crossings that can or 
cannot be resolved depending on the strength of the in- 
teraction. This leads naturally to the distinction between 
the various regimes of interaction strength. The actual 
calculation of the transport is carried out first for a single 
particle LZ crossing in Sections 7, and in Appendix C. 
This is used as a building block for the calculation of the 
stirring in Sections 8,9. The conclusions and perspectives 
are summarized in Section 10. 



II. THE BOSE-HUBBARD TRIMER MODEL 

The simplest model that captures the physics of quan- 
tum stirring is the three-site Bose-Hubbard Hamiltonian 
(BHH) [HGil Ell (see Fig. 1). This minimal model con- 
tains all the generic ingredients of large BHH lattices and 



therefore often is used as a prototype model in many re- 
cent studies 0, [l^, Ell ■ A classical analysis of this model 
has been performed in [l^, EH where it was shown that 
for appropriate system parameters and initial conditions 
chaotic dynamics would emerge. In this work we con- 
sider adiabatic driving of the ground state preparation 
and therefore chaotic motion is not an issue. 

We consider the site index of the quantum trimer tak- 
ing values i = 0, 1,2. The i=0 site has a potential en- 
ergy vo = e and is regarded as a "shuttle" , while the 
i=l,2 sites are regarded as a two level "canal" (with 
vi=V2=0) . The corresponding ./V-boson BHH is: 

2 u 2 

H = Vihi + jYl ~ 1 ) _ k c (b\b 2 + b\b\) 

1 = 2 = 

-ktiblh + b\b Q ) - k 2 {b\b 2 + b%). (1) 

Without loss of generality we use mass units such that 
h=l, and time units such that intra canal hopping ampli- 
tude is k c =l. Accordingly the two single particle levels 
of the canal are e± — ±1. The annihilation and cre- 
ation operators bi and b\ obey the canonical commuta- 
tion relations = <5j.j while the operators hi = b\bi 
count the number of bosons at site i. The interaction 
strength between two atoms in a single site is given by 
U = 4Trh 2 a s V c{ f/m where V B{t is the effective volume, m is 
the atomic mass, and a s is the s-wave scattering length. 

The couplings between the shuttle and the two ends of 
the canal are k\ and k 2 . We assume that both are much 
smaller than k c Their inverse 1 /k\ and l/k 2 are like bar- 
rier heights, and changing them is like switching valves 
on and off. It is convenient to define the two control 
parameters of the pumping as 

Xl = {h-v)' x ->- (2) 

By periodic cycling of the parameters (Xi,X 2 ) we can 
induce a circulating current in the system. We further 
discuss this controlled stirring process in the next section. 



III. STIRRING 

By periodic cycling of the parameters (X\, X 2 ) we can 
imitate a classical peristaltic mechanism and obtain a 
non-zero amount (Q) of transported atoms per cycle. 
During the driving cycle the total number of bosons re- 
mains constant. The energy is not a constant of motion, 
but in the adiabatic limit considered here, the system 
returns to the same state at the end of each cycle. 

The pumping cycle is illustrated in Figs. 1-3. Initially 
all the particles are located in the shuttle which has a 
sufficiently negative on-site potential energy (X 2 < 0). In 
the first half of the cycle the coupling is biased in favor 
of the ki route (X\ > 0) while X 2 is raised until (say) 
the shuttle is empty. In the second half of the cycle the 
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FIG. 2: (Color online) Scheme of the avoided crossings for 
the one particle problem. Left panels: The two-site system 
(|36|) is prepared with one particle in the ground state which 
initially corresponds to having the particle occupying the left 
site. As the potential e is raised the particle encounters an 
avoided crossing. Due to the slowness of the driving it stays 
in the ground state which implies an adiabatic passage to 
the right site. Right panels: The three-site trimer system is 
prepared with one particle in the ground state which initially 
corresponds to having the particle occupying the shuttle site. 
As the potential e is varied the particle encounters avoided 
crossings with the lower canal orbital. A full stirring cycle 
consists of an adiabatic passage through k\ in the first half 
of the cycle, and another adiabatic passage through &2 in the 
second half of the cycle. See the text for further explanations. 



coupling is biased in favor of the hi route (X\ < 0), while 
X2 is lowered until the shuttle is full. Assuming U=0, 
the shuttle is depopulated via the hi route into the lower 
energy level e_ during the first half of the cycle, and re- 
populated via the hi route during the second half of the 
cycle. Accordingly the net effect is to have a non-zero Q. 

If we had a single particle in the system, the net effect 
would be to pump roughly one particle per cycle. If we 
have N non-interacting particles, the result of the same 
cycle is to pump roughly N particles per cycle. We would 
like to know what is the actual result using a proper quan- 
tum mechanical calculation, and furthermore we would 
like to investigate what is the effect of the interatomic 
interaction U on the result. 

The above description of the stirring process might 
look convincing, but in fact it does not hold in the quan- 
tum mechanical reality. The quantum stirring process is 
in general not a peristaltic process but rather a coherent 
transport effect. This point is best clarified by observ- 
ing that the simple minded picture above implies that 
the amount of pumped particles per cycle is at most N. 
This conclusion is wrong. We shall explain in the next 
section that in principle one can get Q 3> N per cycle. 
The proper way to think about the quantum stirring pro- 
cess is as follows: Changing a control parameter, say X2 
with some constant rate X2 , induces a circulating current 
in the system. Each particle can encircle the system more 
than once during a full cycle. Hence Q ^> N is feasible. 
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FIG. 3: Analysis of the pumping cycle for N particles. See 
the text for further details. For a large cycle that encircles 
the whole shaded region we have Q ~ N. The position of the 
"monopoles" is depicted by black dots: (a) no interactions 
(all monopoles are "piled up" at the same position) (b) with 
interactions. In panel (c) we schematically plot the energy 
levels along the Xi = axis for a system corresponding to 
TV = 3 bosons. Note that energy levels that correspond to 
non-participating states (those with non-zero occupation of 
the upper canal orbital) are not plotted. 



IV. THE ADIABATIC PICTURE 

In analogy with Ohm's law (where X is the magnetic 
flux, and —X is the electro motive force), the current 
is / = — G\X\ if we change X\ and / = — G2X2 if we 
change X2, where G\ and G2 are elements of the geo- 
metric conductance matrix. Accordingly 



Q = I Idt=- I(G 1 dX 1 + G 2 dX 2 ). 



(3) 



In order to calculate the geometric conductance we use 
the Kubo formula approach to quantum pumping I17I1 
which is based on the theory of adiabatic processes [18j . 
It turns out that in the strict adiabatic limit G is re- 
lated to the vector field B also known as "two-form" in 
the theory of Berry phase. Namely, using the notations 
Bi = — G2 and B2 = G\ we can rewrite Eq.Q as 



Q = d> B • ds, 



(4) 



where we define the normal vector ds = (dX^, ~dXi) as 
illustrated in Fig. 3. The advantage of this point of view 
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is in the intuition that it gives for the result: Q is related 
to the flux of a field B which is created by "magnetic 
charges" in X space. For U=0 all the magnetic charge 
is concentrated in one point. As the interaction U be- 
comes larger the "magnetic charge" disintegrates into TV 
elementary "monopoles" (see Fig. 3). In practice the cal- 
culation of B is done using the following formula: 



E 

n^riQ 



2 Im[X non ] 

(E n -E no y 



(5) 



where the current operator is conveniently defined as 



h(blbx - b%) + k 2 (b% - blb 2 ) 



while the generalized force operators arc defined as 

dU 



T 3 = - 



dX 4 



(6) 



(7) 



and are associated with the control parameters Xj. The 
index n distinguishes the eigenstates of the many-body 
Hamiltonian. We assume from now on that n is the 
BEC ground state. 



V. THE TWO-ORBITAL APPROXIMATION 

We assume an adiabatic process that involves only two 
orbitals: the shuttle orbital and the orbital of the lower 
canal level (see right panel of Fig. 2). Accordingly we can 
simplify the Hamiltonian in a way that illuminates the 
physics of the quantum stirring process and simplifies the 
formal treatment. For pedagogical reasons we consider 
first the single particle (TV=1) case, and extend the cal- 
culation to the many-body case in the next section. The 
model Hamiltonian in the position basis is: 




H = 



The current can be measured on the 
the 2 i — > bond, accordingly: 



(8) 



1 bond or on 



1: 



2^0 — 



-iki 
iki 


ik 2 

~ik 2 



(9) 
(10) 



For zero couplings (ki = k 2 = 0) the eigenenergies of the 
orbitals are e± = ±1 and £q = e. The corresponding 



eigenstates are 





(11) 
(12) 



The Hamiltonian and the current operators in this orbital 
basis are: 



H = 




(13) 



and 



2qi-»1 = 



-2>->0 



ki + k 2 



k\ + k 2 



- 


-in 


—in 


IK 








in 











IK 


—in 


— IK 








IK 









(14) 
(15) 



where the effective coupling between the shuttle orbital 
and the lower canal orbital is 



fei + k 2 
k = -= — 

and the net "splitting ratio" is [25 

ki - k 2 



X 



ki 



(16) 



(17) 



The stirring cycle starts with all particles localized in the 
shuttle orbital and the adiabatic particle transport takes 
place during the avoided crossing of this orbital with the 
lower canal orbital. Accordingly, we focus on the upper 
left (2 x 2) submatrix: 



n 



S —K 
-K -1 



(18) 



In practice it is more convenient to define the averaged 
current operator X = (2oi->i + 2"2i->o)/2 whose matrix rep- 
resentation is 



1 = 



— IK 

in 



(19) 



The advantage of this definition is that within the two 
halves of a symmetric pumping cycle, while e is raised 
or lowered, the same amount of particles is being trans- 
ported. Thus in order to get Q for a full cycle, we simply 
double the integrated current over a half cycle, where the 
variation of the control parameter e is monotonic start- 
ing from a very negative initial value. For completeness 
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FIG. 4: (Color online) Evolution of the energy levels, the site occupation and the conductance for N = 16 particles and 
fci = 2 x 1CP 4 , = 1 x 10~ 4 . We refer to three representative values of U, which are indicated on top of each set of panels. 
Upper panels: the lowest N + 1 energy levels E n which dominate the conductance Gi are plotted as a function of X2 = s. 
The insets represent magnifications of the indicated areas. Middle panels: the site occupations no (blue A), m (black o), ri2(red 
□ ). Note the steps of size 1 for the dot-occupation and 1/2 in the wire-sites occupation in subfigure (c). Lower panels: the 
corresponding conductance G2 as a function of e. Numerical results are represented by solid black lines while the dotted red 
line corresponds to the analytical result (|47|l in (a) and to (|55p in (b), (c). 



we also write the matrix representation of the generalized 
force that is conjugate to X 2 = s 

^=(;s)- ( 2 °) 

The above expressions for TL and X and T completely 
define the transport problem in the case of one parti- 
cle. Within the two orbital approximation the calcula- 
tion of Q therefore reduces to the study of a single LZ 
crossing, which we discuss in Section VII. 



VI. THE ADIABATIC VARIATION OF THE 
MANY-BODY LEVELS 

We now turn to the analysis of the many-body prob- 
lem. The first step is to understand the evolution of the 
eigenenergies E n as e is varied, while the other param- 
eters are kept constant. It is convenient to rewrite the 



BHH as 

H = H Auta M + W C ana, + Kp. (21) 

H Bhuttlc = — n (n -l) + en 

= ^Mni-l)+n 2 (n 2 -l)]-(i4&i + &I&2) 

while the operators for the current and the generalized 
force are 

2o^i = iki{b\b Q - b\bi) 
l 2 ^o = ik2{blb 2 - b\b ) 

f = -no- (22) 

In what follows we assume < k\,k 2 *C k c = 1 , and 
N\U\ ^ k c = 1, and consequently generalize the "two 
orbital approximation" of the previous section to the case 
of N > 1 particles. 

In the zeroth order approximation k± and k 2 are ne- 
glected; later we take them into account as a perturba- 
tion. For k\ = k 2 = the number (n) of particles in the 
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FIG. 5: (Color online) Evolution of the energy levels, the site occupation and the conductance for strong attractive interaction 
U = — 1 (see Fig. U for legend and parameters). The particles are "glued together" and roll like a classical ball from the shuttle 
to the left wire site as can be seen from the middle panel where the site occupations no(blue A), ni(black o), «2(red □) are 
plotted. The width of the transition is exponentially small in N, as explained in the text, and therefore cannot be resolved 
numerically. 



shuttle becomes a good quantum number. The other 
N—n particles occupy the lower orbital (e_) of the canal 
because we assume NU <C k c . Hence the many-body 
energies are 



E n = E BhMle (n) + E c ^(N-n), 
where n=0, 1, .., N, and 

Puttie = en+^U(n-l)n 



(23) 



(24) 



E... 



-(N—n) + -U(N—n—l)(N—ri). (25) 



From the degeneracy condition E„ — E n -\ = 
(n=l,2, ...N is the number of particles in the shut- 
tle) we determine the location of the n i— > (n— 1) crossing 
to be 



-1 + -Ux (N-3n+2). 



(26) 



Accordingly, we conclude that the N crossings are dis- 
tributed within 

-1-(JV-1)Z7 < e < -1 + -{N-l)U. (27) 

One can introduce a rescaled control variable e which 
reads 

£ + 1 



(N- 1)U' 



(28) 



and its support is —1 < e < 1/2. The distance between 
the crossings, while varying the shuttle potential s, is 
(3/2) U. Once we take n into account we get avoided 
crossings, whose width we will estimate in the next para- 
graph. 

Within the framework of the two-orbital approxima- 
tion, the truncated many-body Hamiltonian matrix takes 
the form 

^nm — E n 8 n ^ ra ^n^n,n±l7 (^^) 

where n=0, ■ ■ ■ , N and the couplings are defined as 
K n = (n—l\Tl\n). For example, in the iV=3 case we have 

/ E -Ki 

— Kl El —K2 
— K2 Ei — K3 
V -K 3 E 3 

The calculation of K n involves the matrix elements of 
6j6oj leading to 

K n = [{N+ 1 -n)n] 1/2 K. (31) 

An analogous expression applies to the truncated current 
operator: 








n = 



(30) 
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(32) 
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FIG. 6: (Color online) Integrated density (IDoS) of avoided 
crossings ("magnetic monopoles") for various values of the 
parameters U , ki , fe and the boson number N as a function 
of the rescaled on-site potential £. The support of the IDoS 
is predicted by Eq. (|28[1 to be £ = [—1,0.5] which is nicely 
confirmed. 



For large U , as e is varied, we encounter a sequence of 
distinct LZ transitions: 



|3) 



|2) 



|1) 



(33) 



The distance between avoided crossings is of order U 
while their width is 



FIG. 7: (Color online) Scaling behavior of the conductance G 
in the gradual crossing regime (U/k « 4.7). We have N = 16 
particles, and the curves correspond to various values of k 
and A. The x-axis is the rescaled control variable £ ()28p and 
the j/-axis is scaled accordingly to preserve the net charge. 
Additionally, the ordinate was scaled by the expected charge 
Q for a half-cycle according to Eq. ((48} and indeed the area 
under the curve is Q « 0.5. One observes that the curves fall 
on top of one another which confirms the dependence of G 
on the ratio U/k. Additionally, we overplotted the theoret- 
ical expression (solid orange line) for the conductance (|55|) . 
Although this expression is not expected to be valid in the 
intermediate regime, the agreement is pretty good. Also the 
agreement with the estimation from Eq. (|56[1 corresponding 
to a constant value of — G ~ 0.315 is apparent. 



Se n = K n . (34) 

The widest crossings are at the center with Se n ~ Nn. 
This width should be smaller than the spacing U between 
avoided crossings, else they merge and we no longer have 
distinct crossings. The other extreme possibility is to 
regard U as the perturbation rather than n. The width of 
the one-particle crossing is k, and it would not be affected 
by the many-body interaction as long as the span NU is 
much smaller than that. We therefore deduce that for 
repulsive interaction there are three distinct regimes: 

U <C k/N mega crossing regime 
k/N < U < Nk gradual crossing regime (35) 
U 3> Nk sequential crossing regime 

Accordingly, depending on the ratio U/k we expect dif- 
ferent results for G2(X). Indeed in a later section this 
expectation is confirmed both analytically and numeri- 
cally (Figs. 0][5]). In Fig. [5] we report the integrated den- 
sity of avoided crossings (IDoS) for various values of U, 
k, A and number of bosons N. We find that all points fall 
in the predicted range confirming nicely the scaling rela- 
tion l|28p. If U/k is small these avoided crossings merge 
and cannot be resolved. In a later section we discuss the 
implied scaling relation for the conductance (Fig. [7]). 



VII. TRANSPORT DURING A LZ CROSSING 

The prototype example for an adiabatic crossing is the 
LZ problem. In this section we discuss the analysis of the 
transport during a LZ crossing, while in the next section 
we shall use the obtained result as a building block for the 
analysis of the transport during a stirring process. The 
following treatment assumes a strict adiabatic process. 
A more advanced treatment that takes into account non- 
adiabatic transitions can be found in p2j , where also the 
resulting fluctuations in Q are calculated. 

Consider a single particle in a two-site system (see left 
panels of Fig. 2). The coupling between the sites is re 
while the on-site potentials are = ±e/2, where i = 1,2 
labels the sites. The Hamiltonian is 

n, = ( £ H _-; /2 ) . (36) 

At time t = — oo, the control parameter X = e is very 
negative and the particle is on the left site (i = 1). 
Then X is increased adiabatically and the levels experi- 
ence an avoided crossing leading to the adiabatic transfer 
of the particle to the right site (i = 2). The current op- 
erator and the generalized force operator are represented 
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by the matrices 



and 



— in 
Mj ~ I in 



T--\ 1/2 ° 
tJ ~ \ -1/2 



(37) 



(38) 



The instantaneous eigenenergies of the LZ Hamiltonian 
are labeled as n = — (lower) and n=+ (upper): 



and the associated cigenstates are 

1^-) = 
\E + ) = 

where 



(39) 



cos(6>/2) A 
sin(0/2) ) 

-sin(0/2) A 
cos(0/2) J 



n = Ve 2 + (2k) 2 
9 = — arctan (2n/e) 



(40) 



(41) 
(42) 



Note that 8=0 at t= — oo evolves to 6=ir at t=oo. The 
matrix representation of X and T in this basis is 



Xn 



— IK 
IK I ' 



and 



J n 



1 / cos(0) -sin(0) 



2 V - sin(0) - cos(0) 



(43) 



(44) 



Now we can use Eq.© to obtain the geometric conduc- 
tance: 



G(s) 



2k 2 



ft 2 [e 2 + (2k) 2 ] 3 / 2 ■ 



(45) 



Since we assume here a strictly adiabatic process we ex- 
pect 100% transfer efficiency, and indeed: 



Q 



[I)dt = - Gds 



= 1. 



(46) 



For pedagogical reasons, we rederive the result (j43]) for 
the conductance G(e) in a straightforward way from the 
Schrodingcr equation in Appendix C. This might add to 
the understanding of the general discussion of quantum 
stirring presented in a later section. 



VIII. TRANSPORT DURING STIRRING: 
THE CASE U = 

We have realized that the stirring problem of one par- 
ticle reduces to the LZ-like crossing problem Eqs. (|36l37[) . 
Disregarding the different on-site energies, the Hamilto- 
nian is the same as in the LZ-problem Eq. (|36|) . The sig- 
nificant difference is related to the current operator. Its 
multiplication by the splitting ratio A is the fingerprint of 
the non-trivial topology. We further discuss the physics 
behind A at the end of this section. On the technical side, 
all we have to do is to multiply the result obtained in the 
previous section by this factor. If we have N non inter- 
acting Bosons, the result should be further multiplied by 
N, leading to 



g{x 1 ,x 2 ) = -\n\- ^ 

2 [( £ _£_)2 + 



(2k) 2 ] 3 / 2 ■ 



(47) 



In this expression e is determined by X2 and A is deter- 
mined by X\ , while k is conveniently regarded as a fixed 
parameter. For a full cycle we get 



Q = f(I)dt = NX. 



(48) 



We note that this result assumes a symmetric stirring 



cycle such that ki 



and &2 = fcsmaii in the first half 



of the cycle, while k\ = fc sman and k 2 = A: largo in the second 
half of the cycle. Accordingly 



= N 



(49) 



It is more illuminating to denote the Jfi-radius of the 
pumping cycle by R , as illustrated in Fig. [3] and to re- 
write the latter expression as follows: 



Q = N 



[1 + {kR) 2 ] 1 ' 2 - 1 



kR 



(50) 



In particular for small cycles we get a linear dependence 
on the radius 

Q w NkR, (51) 
while for large cycles we obtain the limiting value 



Q 



N. 



(52) 



In a two-site topology the amount of particles that are 
transported during a strictly adiabatic LZ crossing is ex- 
actly N. In contrast to that, in the stirring problem that 
we study, we have non trivial "ring" topology, and there- 
fore Q is multiplied by the splitting ratio A. In practice 
ki and k 2 are positive and accordingly |A| < 1. But in 
principle we can have |A| > 1. This would happen if ki 
and fc 2 were negative. In such case the lower orbital of the 
"canal" is anti-symmetric rather than symmetric. Going 
through the derivation one observes that the same results 
apply with 



A 



1/A. 



(53) 



For small cycles we find 

Q ps N[kR}-\ (54) 

which means that we can circulate Q N particles per 
cycle. This demonstrates our statement that quantum 
stirring is not a classical-like peristaltic process, but 
rather a coherent transport effect. 

The results above become more transparent if we no- 
tice that they reflect the topology of the (Xi,X 2 ,X 3 ) 
space, where X 3 is a fictitious Aharonov-Bohm flux to 
which the operator 2 is conjugate [see Appendix A]. In 
this extended space the field B that appears in Eq.Q 
has zero divergence, with the exception of the "Dirac 
monopolcs" which are located at points where E no has a 
degeneracy with a nearby level. If the shuttle orbital and 
the lower orbital of the "canal" have the opposite parity 
this degeneracy point is located at X = (0, £_,0) which 
is in the plane of the pumping cycle. But if the shut- 
tle orbital and the lower orbital of the "canal" have the 
same parity, then this degeneracy point is displaced off- 
plane. Accordingly we get the divergent result Eq.([51|) 
or the non-divergent result Eq.(|5ip. It is important to 
realize that because of the gauge invariance under the 
transformation X 3 t— > A3+27T the degeneracy point is 
duplicated, leading to a "Dirac chain". In the far field 
this chain looks like a charged line, and accordingly in 
the far field we get Eq. (|5"12|) which in leading order is not 
sensitive to the parity of the orbitals. 

If we had only one particle, the degeneracy at the cen- 
ter of the pumping cycle would correspond to a single 
Dirac monopole. If we have N non-interacting particles 
we have in fact N Dirac monopolcs at the same location. 
As the interaction U is turned on this "pile" disintegrates 
into N elementary "monopoles" (see Fig. 3). This will 
be further discussed in the next section. 



IX. TRANSPORT DURING STIRRING: 
THE CASE U 

If we have a rectangular pumping cycle in A-space, it 
can be closed at X 2 = ±oo, where the influence of the 
change in X\ can be safely neglected since the monopolcs 
are located on the X\ = axis around e ps — 1. Accord- 
ingly, the predominant contribution to Q results from 
the dX 2 variation and therefore we refer from now on 
to G(e) = G2(X) only. An overview of the numerical re- 
sults for the conductance is shown in Figs. (OO, where 
we plot G as a function of X 2 = e for various interaction 
strengths U . Besides G we also plot the ^-dependence 
of the energy levels and of the site population. Five rep- 
resentative values of U are considered including also the 
case of weak/strong attractive interactions U < 0. Fi- 
nally, in Fig. [5] we plot the implied dependence of Q on 
the radius of the pumping cycle, again for several repre- 
sentative values of U. 




_l i I i I i L_ 

0.8 0.9 1 1.1 
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FIG. 8: (Color online) Dependence of Q on the X2-radius R 
of the pumping cycle for the system presented in Fig. [4] 
The curves correspond to different values of the interaction 
strength U and represent the net amount of particles (Q) 
transported during the first half of a symmetric pumping cy- 
cle centered around X2 = e = 0. For vanishing interactions 
U (dotted line) all the particles are transported once the cycle 
encloses the value X2 = — 1 (mega-crossing). As U is being in- 
creased the transport gradually changes leading eventually to 
a step-like behavior (solid line) for large repulsive interactions 
(sequential crossing). 



A. Mega crossing 

For small positive values of U, the dynamics appears 
to be the same as in the U — case. Namely, all the 
particle cross "together" from the shuttle orbital to the 
E- canal orbital. We call this type of dynamics "mega 
crossing". In the language of the adiabatic picture this 
means that all the Dirac monopoles are piled up in the 
center of the pumping cycle. Once the interaction U 
is turned on this "pile" disintegrates into N elemen- 
tary "monopolcs" . Consequently the dependence of Q on 
the radius of the pumping cycle becomes of importance. 
Disregarding fluctuations that reflect discreteness of the 
magnetic charge, the value of Q is determined by the 
number of monopoles that are encircled by the pumping 
cycle. Thus, by measuring Q versus R we obtain infor- 
mation on the distribution of the monopoles and hence 
on the strength of the interatomic interactions. 



B. Sequential crossing 

In the sequential crossing regime the transport can be 
regarded as a sequence of LZ crossings. Each LZ crossing 
is formally the same as the LZ crossing of a single particle 
in a two-site system, while the topology is reflected by 
the splitting ratio A and the interaction is reflected in the 



10 



scaled coupling constants K n . Accordingly we get 

N 

2 A 



G 



1 N 



2ni 



J [(£-£„) 2 + (2 K „)2]3/2' 



(55) 



We overplot this formula in the lower panel of Fig. 5J; 
where an excellent agreement is observed. 



C. Gradual crossing 

For intermediate values of U (weak repulsive interac- 
tion), i.e. in the range k/N <C U <C Nk, we find neither 
the sequential crossing of Eq. (|55p , nor the mega-crossing 
of Eq. ([47]). but rather a gradual crossing. Namely, in 
this regime, over a range AX2 = (3/2)(N—l)U we get a 
constant geometric conductance: 



G 



-X- 



1 

'SU 



(56) 



which reflects in a simple way the strength of the inter- 
action. This formula has been deduced by extrapolating 
Eq. (|55[) . and then was validated numerically (see lower 
panel of Fig. SJs). In this regime it is more illuminating 
to plot the scaled conductance 



G(i) 



(N-1)U 



G(s) 



(57) 



which is implied by the scaling (|28[) of e. The numerical 
results are reported in Fig. [7] The shape of the plot 
depends only on the dimcnsionless parameters U /k and 
N. The curves corresponding to different A and U values 
(but with the same constant ratio U/n) fall nicely one 
onto the other with good accuracy, confirming the scaling 
relation (T5"7) . 



D. Attractive interaction 

So far we have discussed repulsive interactions for 
which the iV-fold "degeneracy" of the U = mega cross- 
ing is lifted and we get a sequence of N avoided crossings. 
Also for U < this TV-fold "degeneracy" is lifted, but in 
a different way: The levels separate in the "vertical" (en- 
ergy) direction (see upper panels of Fig. [5]) rather than 
"horizontally" (see upper panels of Fig. [4]). 

In the U < regime all the particles execute a sin- 
gle two-level transition from the shuttle to the canal 
(see Fig. [S£i). This transition happens directly from the 
n = TV state to the n = state. The coupling between 
these two states is exponentially small in TV because it 
requires a virtual TV-th order transition in perturbation 
theory via the intermediate n values. 

For sufficiently strong attractive interactions 
(|TV[/| 3> 1) all the particles arc glued together and 
behave like a classical ball that rolls from the shuttle to 
one of the canal sites (see Fig. [SJd). When the sign of 



Xi is reversed the ball rolls from one end of the canal 
to the other end (not shown). This has to be clearly 
distinguished from the TV-fold degenerated transition to 
the lower canal level which is observed in the [7=0 case. 



X. CONCLUSIONS 

In this paper we have considered the transport induced 
in a few site system that contains condensed particles. 
One new aspect that has emerged throughout the analy- 
sis is the existence of four distinct dynamical regimes de- 
pending on the strength of the interatomic interactions. 
This observation also applies to studies of LZ crossings in 
dimer systems and has further implications regarding the 
quantum stirring in topologically non-trivial systems. It 
should be clear that the analysis of our "driven vortex" 
requires the toolbox of adiabatic processes, and it should 
be distinguished from the ignited stirring of Refs. [23l.[2^|. 

The actual measurement of induced neutral currents 
poses a challenge to experimentalists. In fact, there is 
a variety of techniques that have been proposed for this 
purpose. For example one can exploit the Doppler effect 
at the perpendicular direction, which is known as the 
rotational frequency shift [25| . 

The analysis of the prototype trimer system reveals the 
crucial importance of interactions. The interactions are 
not merely a perturbation but determine the nature of 
the transport process. We expect the induced circulating 
atomic current to be extremely accurate, which would 
open the way to various applications, either as a new 
metrological standard, or as a component of a new type 
of quantum information or processing device. 
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APPENDIX A: THE BERRY PHASE AND 
THE B-FIELD 

In this Appendix we explain how the field B emerges in 
the theory of Berry phase. In the presentation below we 
follow the notations as in [l7| . Given a time dependent 
Hamiltonian TC(X(t)) where X = (Xi, X2, X3) are con- 
trol parameters it is convenient to expand the evolving 
wavefunction in the adiabatic basis 



*(*))= J2 a n(t) "(*(*)) 



(Al) 



Then the Schrodinger equation id t \t^) = H(X(t)) be- 
comes 

da n 



iE n a n + i ^2 ^2 XjA 3 nm a ri 



(A2) 
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where we defined 



(A3) 



Differentiation by parts of dj(m(X)\n(X)) = implies 
that A 3 nm is a Hermitian matrix. We denote its (real) 
diagonal elements as 



A\X) ee Ai, 



(A4) 



The line integral over the vector field A(X ) along a closed 
driving cycle gives the Berry phase 



BerryPhasc = i> A • dX. 



(A5) 



we find that the off- 



Using d X] (m(X)\H\n(X)) 
diagonal elements of A J nm can be written as 



A 3 



Em E n 



on 



dXi 



, = ( A6) 



The "1-form" A? is formally like a vector potential, and 
we can associate with it a gauge invariant "2-form" B k i 
which is formally like a magnetic field: 



d k A> 



d 3 A k 



—2lm(dkn\djn) 
— 21m ^ ' A nm A J mn 



E 



2Im[^ m ^U 

(E m -E n y 



(A7) 



In order to make the magnetic field analogy more trans- 
parent we assume that we have three control parame- 
ters {Xi, X%, A3). Then it is natural to associate with 
the antisymmetric matrix B 3 a field whose components 
are Bi = —B 32 and B 2 = B 31 and B 3 = B 12 . This 
B(A) field has has zero divergence everywhere with the 
exception of the Dirac monopolcs at points of degeneracy 
Dirac monopoles have quantized charge such that their 
flux is an integer multiple of 2tt. The Dirac monopoles 
must be quantized like that, otherwise Stokes' theorem 
would imply that the Berry phase is ill-defined. 



APPENDIX B: THE KUBO FORMULA AND 
THE B-FIELD 

The Kubo formula is traditionally used in order to 
calculate the response of a driven system in the linear 
response regime. Given a time dependent Hamiltonian 
Tt(X(t)) the linear DC-response of the system is ex- 
pressed as 



(T k 



(Bl) 



where the generalized conductance matrix G kj can be 
calculated from the Kubo formula. This matrix can 
be decomposed in a symmetric and an anti-symmetric 
part which account for the dissipative and non-dissipativc 
effect of the driving respectively. Here we consider a 
strictly adiabatic driving: Although the energy is not a 
constant of motion the system returns to the initial state 
at the end of each cycle. In this case there is no dissipa- 
tion and thus we consider only the anti-symmetric ("ge- 
ometric" ) part of G k] which is identified as B kj . We fur- 
ther illuminate this identification in the next appendix. 

In the stirring problem there are two control parame- 
ters which we call X\ and A 2 , and the current operator 
X = J- 3 is conveniently regarded as conjugate to a ficti- 
tious Aharonov-Bohm flux parameter A3. If the particles 
were charged we could regard A3 as an actual control pa- 
rameter; then A3 would be the electro-motive-force and 
G 33 would be the conventional Ohmic conductance. In 
Section 4 we use simplified indexing, namely Gj = G 3j . 
Accordingly Eq. (|B1[) leads to Eq. ^ and in the adiabatic 
limit we obtain Eq.Q with Eq.© which follows from 
Eq.flAT]). 



APPENDIX C: ALTERNATIVE DERIVATION 
FOR THE GEOMETRIC CONDUCTANCE 

The definition of B in Appendix A illuminates its ge- 
ometric interpretation in the context of the the Berry 
phase formalism, but does not help in understanding why 
it emerges in the linear response calculation as described 
in Appendix B. For this reason it might be helpful to de- 
rive Eq. (|4"5)) directly from the Schrodinger equation. For 
the LZ crossing problem it becomes 



da n _ _. p 
at 



nm^m j 



(Cl) 



where e is a small parameter. Using the explicit expres- 
sions Eq. (|4"0|) for the adiabatic eigenstates, the definition 
Eq. (|A3j) . and the identities d x arctan(a;) = l/(l+a; 2 ) and 
39/ de = 2n/n 2 we find 



A Tl 



K 



-i 

1 



(C2) 



The zero order adiabatic eigenstates of the Hamiltonian 
H = H - eA are |£L) and \E+) of Eq.gD}. The first 
order adiabatic eigenstates are found from perturbation 
theory. In particular the lower state is 



= \E_ 



(C3) 



It is important to realize that the expectation value of 
the current operator vanishes in zcroth order. This is be- 
cause the the zeroth order Hamiltonian Tl is time-reversal 
symmetric. It is the i perturbation term that breaks the 
time-reversal symmetry leading to 



{-\T\-) = 2e 



n 3 ' 



(C4) 
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Using the notation (J) = — Ge we see that this result is 
in agreement with Eq. (|45[) as expected. 
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